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ABSTRACT 

We present observations of the afterglows and host galaxies of three short-duration gamma-ray bursts (GRBs): 
100625 A, 101219A and 1 101 12A. We find that GRB 100625A occurred in a z = 0.452 early-type galaxy with a 
stellar mass of w 4.6 X IO^Mq and a stellar population age of « 0.7 Gyr, andGRB 101219A originated in a star- 
forming galaxy at z = 0.718 with a stellar mass of w 1.4 x 10^ Mq, a star formation rate of « 16 Mq yr"', and 
a stellar population age of w 50 Myr We also report the discovery of the optical afterglow of GRB 1101 12 A, 
which lacks a coincident host galaxy to / > 26 mag and we cannot conclusively identify any field galaxy as 
a possible host. From afterglow modeling, the bursts have inferred circumburst densities of « 10"^- 1 cm"-', 
and isotropic-equivalent gamma-ray and kinetic energies of w 10^*'- 10^' erg. These three events highlight the 
diversity of galactic environments that host short GRBs. To quantify this diversity, we use the sample of 36 
Swift short GRBs with robust associations to an environment 1/2 of 68 short bursts detected by Swift to May 
2012) and classify bursts originating from four types of environments: late-type (« 50%), early-type (w 15%), 
inconclusive (w 20%), and "host-less" (lacking a coincident host galaxy to limits of > 26 mag; « 15%). To 
find likely ranges for the true late- and early-type fractions, we assign each of the host-less bursts to either 
the late- or early-type category using probabilistic arguments, and consider the scenario that all hosts in the 
inconclusive category are early-type galaxies to set an upper bound on the early-type fraction. We calculate 
most likely ranges for the late- and early-type fractions of w 60 - 80% and « 20 - 40%, respectively. We 
find no clear trend between gamma-ray duration and host type. We also find no change to the fractions when 
excluding events recently claimed as possible contaminants from the long GRB/collapsar population. Our 
reported demographics are consistent with a short GRB rate driven by both stellar mass and star formation. 
Keywords: gamma rays: bursts 



1. nvlTRODUCTION 

Observations of the galactic environments of cosmic 
explosions provide invaluable insight into their underly- 
ing progenitor populations. For example. Type la su- 
pernovae (SNe) originate in both star -forming and ellip- 
' tical galaxies dO emler & Tinslev 1979; _[yan den Bergh et al.l 

■ l2005l;lMannucci et al.ii2005t:iLi et al.ii201 ll) ' consistent with an 
] evolved progenitor and a n event rate that trac es both stellar 

mass and star formation (ISullivan et al.ll2006h . In contrast, 
SNe of types II and Ib/c are found to occur only in spiral and 
irregular galaxies, indicating that these e vents result from the 
core-collapse of young, massive stars ("van den Bergh et all 

■ |1005; Hakobyan et al. 2008; Li et al. 201 1) a nd a rate tracing 
recent star formation (iKellv & Kirshneil2012t I Anderson et al.l 
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In the c ase of long-durati on gamma-ray bursts, (GRBs; 
T()o > 2 s; iKouveliotou et alj ri993) the link to star-forming 
host galaxies helped to establish that their progenitors are 
massive stars CDiorgovski et al. 19 98; Le Floc'h et al. 20031 
iFruchter et al.l 12006: Wain wright et al. 2007a). Furthermore, 
a decade of concerted efforts to characterize the stellar pop- 
ulations of long GRB hosts revealed young stellar popu- 
lation ages of < 0.2 Gyr, a mean stellar mass of « 2 x 
1O^M0, and inferred UV/optical star formation rates (SFR) 
of Ki l-50Mf;i vr~' JChri stensen et al. 2004; Sa vaglio et all 
120091; lLeibler&Bergeril2010; .Laskai'et aL2011i) . In addi- 
tion, the spatial locations of long GRBs with respect to their 
host galaxy centers (with a mean of ^ 1 half-light radius; 
iBloom et al. 2002) and their concentration in bright UV re- 
gions of their hosts (iFruchter et all 12006 ) provided a direct 
association between long GRBs and star formation. 

In contrast, the origin of short GRBs (Tgo < 2 s) is less 
clear, as the first few afterglow discoveries led to associations 
with both elliptical (Berger et al. 2005; Castro-Tirado et aH 
12005^, iGehrelset al] [2005; Hiorth et al. 2005a; Bloom et all 
| 2006i) and star-for ming (Fox et al. 2005; Hiorth et al. 2005 W 
ISoderberg et al.l2 006: Grupe et al. 2006; Burrows et al. 20061 
host galaxies, demonstrating that at least some short GRBs 
originate from older stellar populations. Studies primarily fo- 
cused on the sample of bursts with sub-arcsecond localization 
have shown the population of hosts to be dominated by late- 
type galaxies, albeit with lower specific SFRs, higher lumi- 
nosities, and h igher rnetalli cities than the star-forming hosts 
of long GRBs ( lBerged2009l) . Modeling of the spectral energy 
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distributions of short GRB host galaxies has led to a broad 
range of inferred ages, r « 0.03-4.4 Gyr, and an average 
stellar mass of « 2 x lO'^M© jLeibler & BergerilMT O). A 
detailed analysis of their sub-galactic environments through 
Hubble Space Telescope observations has demonstrated that 
on a verage, short GR Bs have offsets from their hosts of « 5 
kpc (iFong et al.ll2010l) . while a growing subse t which lack co - 
incident hosts may have offsets of > 30 kpc (lBergeril2010ab . 
Finally, an examination of short GRB locations with respect to 
their host light distributions revealed that they under-represent 
their host UV/optical light (Fong et al. 2010). These results 
are consistent with theoretical expectations for NS-N S/NS- 
BH mergers (E ichler et aLlll989t iNarayan et al.lll992l) . with 
potential minor contribution from other proposed progeni- 
tors, such as the accretion-induced-collapse of a WD or NS 
dOinet all [1998; Levanetal. 2006b; Metzgeretal. 2008) or 
magnetar flares (iLevan et al.l l2006b; Cha pman et al.ll200 8). 

However, the majority of short GRB host galaxy studies 
published thus far primarily concentrate on bursts with sub- 
arcsecond localization from optical afterglows. While these 
events have the most unambiguous associations with host 
galaxies, the fraction is only ~ 1/3 (23/68 to May 2 012) of 
all sh ort GRBs detected by the Swift satellite ( Gehrels et al.l 
12004ft . The faintness of their op tical af terglows (« 23 mag at 
^ 10 hr after the burst; .Bergerll2010al is likely attributed to 
a combination of a low energy scale (iPanaitescu et al.ll200l1) 
and circumburst densities. Therefore, if there exist correla- 
tions between these basic properties and host galaxy type, 
the selection by optical afterglows may affect the relative 
rates of short GRBs detected in early- and late-type host 
galaxies. An alternative route to sub-arcsecond localization 
is through the X-ray detection of an afterglo w, which does 
not ne cessarily depend on circumburst density (iGranot & Saril 
|2002) with Chandr a; however, only tw o such cases have been 
reported thus far (iFong et al.l 120 12[ iM argutti et al. 2012al 
ISakamotoetal.1120121) 

Demographics which accurately represent the bulk of 
the short GRB population are imperative in understand- 
ing the link to the progenitors. In particular, the late- 
to-early-type host galaxy ratio will infor m whether stellar 
mass or SFR drives the short GRB rate (iLeibler & Bergerl 
120101). and will help to constrain the delay time distribu- 
tion (IZheng & Ramirez-Ruizll2007l) . Furthermore, a recent 
study based on 7-ray properties (spectral hardness and du- 
ration) claims that there is a non-negligible fraction of con- 
taminants from collapsars in the Swift short GRB population 
(|Bromberg et al. 2012). Thus, an examination of how this 
fraction affects the environment demographics will aid in as- 
sessing the true contamination. 

Fortunate ly, the detection of X-ray afterglow s with 
Swift/XRT dGehrels et al.ll2004l: iBurrows et al.ll2005h enables 
positions with ~few arcsecond precision in ss 60% (40/68) 
of all Swift short GRBs. In the majority of such cases, these 
XRT positions coupled with dedicated optical/NIR searches 
for host galaxies haveprovided meaningful associations to a 
galactic environmen{3 While such bursts with XR T posi- 
tions have been studied as single events (e .g . Gehrel set al.l 
|2005; Bloom et al."2006, 2007|lPerley.etaJj|2012), the entire 
sample has not been studied in detail alongside bursts with 
sub-arcsecond locaUzation. 

The large majority of the remaining fx 40% of Swift short GRBs lack 
afterglow follow-up due to observing constraints unrelated to the burst prop- 
eilies; see §5. 



To this end, we present here X-ray and optical/NIR 
observations of the afterglows and environments of three 
short GRB^I3 localized by Swift fXRT, which highUght the 
diversity of their galactic environments: GRBs 100625A, 
1012 19A, and 1101 12 A. We also present the discovery of the 
optical afterglow of GRB 1101 12A. While GRBs 100625 A 
and 1012 19A have robust associations with host galaxies, 
GRB 1 101 12A lacks a coincident host to deep optical limits. 
We describe the X-ray, optical and NIR observations for these 
three events (§2), present their energy scales and circumburst 
densities inferred from afterglow modeling (§3), and host 
galaxy stellar population ages, masses and SFRs extracted 
from spectroscopy and broad-band SEDs (§4). We discuss 
the stellar population characteristics of these three host galax- 
ies compared to previous short GRB hosts (§5). Putting these 
bursts into the context, we undertake the first comprehensive 
study of host galaxy demographics of both sub-arcsecond lo- 
calized and XRT-localized bursts, by investigating the late- 
and early-type host galaxy fractions for the bulk of the short 
GRB population, and compare host galaxy type to 7-ray prop- 
erties (§6). 

Unless otherwise noted, all magnitudes are in the AB sys- 
tem and are corr ected for Galactic extinction in the direc- 
tion o f the burst (iSchlegel etal.lll998l; ISchlafly & Finkbeined 
I2OI ll) . and uncertainties correspond to Ict confidence. We 
employ a standard ACDM cosmology with IIm = 0.27, J7a = 
0.73, and i/o = 71 km s"' Mpc'^. 

2. OBSERVATIONS 

2.1. GRB 100625 A 

GRB 100625 A was detected by three 7-ray satellites on 
2010 June 25.773 UT: th e Burst Alert Telescope (BAT) on- 
board flie Swift sateflite (Gehrels et aP l2004t iHoUand et al.l 
I2010al) . Konus-Wind (Golenetski i et al.l l2010al) an d the 
Gamma-Ray Burst Monitor (GBM) on-board Fermi jBhatl 
(2010). BAT localized the burst to a ground-calculated posi- 
tion of RA=01''03™11.P, Dec=-39°05'29" (J 2000) with an 
uncertainty of 1.0' radius (90% containment; [Holland et al.l 
Ho 10b), and the burst consisted of two pulses with a to- 
tal duration of Tgo = 0.33 ± 0.03s (15-350 keV) and a flu- 
ence of A = (2.3 ± 0.2) x 10"^ erg cm'^ (15-150 keV; 
iHolland et alj|2010bh . FermZ/GBM observations determined 
£peak = 509!^^ k eV and U = (1.32 ±0.05) x 10"'' erg cm'^ 
(8 - 1000 keV; iBhati 1201 Oh . while Konus-Wind observa- 
tions determined Epeak = 418:!:7g^ ^eV and f^ = (8.3 ± 1.5) x 
10"^ erg cm-2 (20-2000 keV; iGolenetskii eTaLl l2010al) . 
Based on the short duration and high iSpeak, GRB 100625A 
can be classified as a short, hard burst. The 7-ray properties 
are listed in Table [T] 

2.1.1. X-ray Observations 

The X-ray Telescope (XRT) on-board Swift began observ- 
ing the field at 6t = 43 s (St is the time after the BAT trig- 
ger) and detected a fading, uncatalogued X-ray source at 
RA=01''03'"10.9P and Dec=-39°05'18.4 " wifli a final po- 
sitional accuracy of 1.8" radius (90% ; iGoad et al.l [20071; 
lEvans et al.ll2009l: iHoUand et al.ll2010bl; TablelB- 

" We present observations of two additional short GRBs, 100628A and 
100702A, both with published SwiftfXKI localizations (see Appendix). We 
show that the XRT afterglow of GRB 100628A is of low significance, while 
the XRT position of GRB 100702A is contaminated, preventing an unam- 
biguous association with a host galaxy. 
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Table 1 

Short GRB Properties 



GRB 


R.A. 


Decl 


Unceit. 




T90 (15-350 keV) 


f-y (15-150keV) 


References 




(J2000) 


(J2000) 


(") 




(s) 


(erg cm"^) 




GRB 100625 A 


oi^os^io.gp 


-39°05'18.4" 


1.8 


0.452 


0.33 ±0.03 


(2.3 ± 0.2) X 10-' 


1 


GRB 101219A 


04''58"'20.49' 


-02°32'23.0" 


1.7 


0.718 


0.6 ±0.2 


(4.6 ± 0.3) X lO-'' 


2 


GRB 1101 12A 


2l''59"43.85^ 


+26°27'23.9" 


0.14 




0.5 ±0.1 


(3.0 ±0.9) X 10-* 


3, This work 



Note. — References: (D IHoUand et aU2010H (2') IKrimm et aU2010al : (3) IBarthelmv et aU20Ill 




10 10 10 

Time after Burst 6 1 (s) 



10 



Figure 1. SwiftlXBJ light curve of GRB 100625A. The triangle is a 3o- upper 
limit. The entire light curve is best fit with a power law characterized by 
ax = -1.45 ± 0.08 (grey dashed line). 

We analyze the XRT data using HEASOFT (v.6. 1 1) and rel- 
evant calibration files. We apply standard filtering and screen- 
ing criteria, and generate a count ra t e light curv e following the 
prescrip tions from Margutti et al.l (120101) and ^Margu tti et alJ 
JloilH). Our re-binning scheme ensures a minimum signal- 
to-noise ratio of S/N = 4 for each temporal bin. To extract 
a spectrum, we use Cash statistics and fit the XRT data with 
an absorbed power law model (tbabs x ztbabs x pow within 
the XSPEC routine) characterized by photon index, F, and 
intrinsic neutral hydrogen absorption column density, A^H.int, 
in excess of the Galactic column density in the direction of 
the burs t, A^h.mw = 2.1 x lO '^ '^cm"^ (typical unc ertainty of 
- 10%; iKalberla et al] 120051: iWakker et alJ l20Tlh . We uti- 
lize the entire PC data set {5t = 60-10^ s), where there is 
no evidence for spectral evolution. Our best-fit spectrum (C- 
stati. = 0.92 for 95 d.o.f.) is characterized by F = 2.5 ±0.2 
and A^H.int < 1 .7 x lO^' cm'^ (3cr) at z = 0.452 (see SUfor the 
redshift determination). Our best-fit parameters are consistent 
with the automatic spectrum fit produced by Page & Holland 
(llOlO). Applying these parameters to the data, we calculate 
the count rate-to-flux conversion factors, and hence the unab- 
sorbed fluxes (Figure[T). 

To quantify the decay rate, we utilize x^-minimization to fit 
a power law to the data in the form Fx{t) oc f"'', with ax as 
the free parameter The entire XRT light curve ((5f w 80- 10^ 
s, PC mode) is best fit with a single power law with index 
ax = -1 -45 ± 0.08 (x?, = 2. 1 for 7 d.o.f.; Figure[D. 

2.1.2. Optical/NIR Observations and Afterglow Limits 

The UV-Optical Telescope (UVOT) on-board Swift com- 
menced observations at (5f = 56 s but no corresponding source 
was found within the XRT position. The 3a limit over 6t ~ 



87-1.2 X 10^ s in the white filter, which transmits over A = 
1600-7000 A (Poole et_ aL 2008), is > 22.6 mag (not corrected 
for Galactic extinction; (Holland et a l. 2010b). Rapid ground- 
based follow-up in the optical and NIR provided early limits 
on the afterglow of / > 22.8 mag at 5t k. 17 min (Suzuki et 
Hojol) and J > 19.4 mag at St « 8.6 hr (iNaitoet al.ii20ldlK 
GROND observations at 5t ~ 12.2 hr place limits of ^ > 23.6 

. 1: 

mag, and riz > 23 mag (INicuesa Guelbenzu et al. 2012). 

We obtained optical observations of GRB 100625 A with the 
Gemini Multi-Object Spectrograph (GMOS) mounted on the 
Gemini-South 8-m telescope, starting at 6t = 12.4 hr in the 
riz filters in poor seeing conditions (Table |2). We analyze 
the data using the IRAF gemini package, and detect a sin- 
gle source within the enhanced XRT error circle in all three 
filters. To assess any potential fading of the source, we ob- 
tained a second set of observations at 6t w 2.6 d, where the 
source is clearly extended. Digital ima ge subtraction using 
the ISIS software package (I AlardI 12000b shows no residuals 
in all three filters (Figure |2}. We therefore place 3 a limits of 
r > 22.6 mag, / > 22.7 mag and z > 22.8 mag on the optical 
afterglow at St « 12.7 hr (Table|2]l. The GMOS zeropoints are 
determined by sources in common with late-time IMACS ob- 
servations (see below), which are calibrated to a standard star 
field at a similar airmass. Our limits match the GROND l imits 
reported at 6t ?s 12.2 hr (Nicue sa Guelbenzu et alj|2012h . 

In addition, we obtained two epochs of 7-band observa- 
tions with the Persson's Auxilliary Nasmyth Infrared Camera 
(PANIC) mounted on the 6.5-m Magellan/Baade telescope at 
St 1.6 and 6.6 d. We analyze the data using standard proce- 
dures in IRAF. Digital image subtraction shows no evidence 
for fading, with a 3a limit of 7 > 23.9 mag (photometrically 
tied to the 2MASS catalog and converted to the AB system) 
at(5f « 1.6d(Table|21l. 

We obtained late-time griz observations of the the field of 
GRB 100625 A with the Inamori Magellan Areal Camera and 
Spectrograph (IMACS) mounted on Magellan/Baade starting 
on 2010 November 14.11 UT. We also obtained /iTj-band ob- 
servations with the FourStar Infrared Camera mounted on 
Magellan/Baade on 201 1 December 07.16 UT (Table|2l). The 
griz zeropoints are calculated using a standard star field at 
a similar airmass, while the /Tj-band zeropoint is determined 
from point sources in common with 2MASS. Our afterglow 
limit and host galaxy photometry are summarized in Table |2] 

We obtained a spectrum of the putative host galaxy with 
the Low Dispersion Survey Spectrograph 3 (LDSS3) mounted 
on the 6.5-m Magellan/Clay telescope on 2011 October 
21.27 UT. A dithered pair of 2700 s exposures was obtained 
with the VPH-ALL grating, which has a wavelength coverage 
of 4000- 10000 A and a spectral resolution of ss 8 A. We used 
standard tasks in IRAF for data reduction, HeNeAr arc lamps 
for wavelength calibration, and observations of the smooth- 
spectrum standard star EGO 1 for flux calibration. We discuss 
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'^B 100625 A- 
Gemini-S/GMOS i-band 
2010 Jun 26.-30 UT^ '- - 5?^ 




2010 Jun 28.40 UT 



G1 






5 arcsec 








28.8 kpc 









Figure!. Gemini-South/GMOS i-band observations of GRB 100625A. The XRT eiTor circle has a radius of 1.8" (90% containment; black). Images are 
smoothed with a 2-pixel Gaussian. Left: St = 0.53 d in poor seeing conditions (fipwHM = 1.9") with a faint host detection. Center: St = 2.63 d with 0.9" seeing. 
Right: Digital image subtraction of the two epochs reveals no afterglow to a 3(j limit of > 22.7 mag. The host galaxy is marked as Gl. 
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Figures. X-ray afterglow light curve of GRB 101219A, including 
5n7/f/XRT observations (red points) and a Chandra/ ACIS-S observation 
(blue point). Triangles denote 3a upper limits. The data over St «i 200- lO'* s 
are best fit with a power law characterized by ax = —1 37 it 0. 1 3 (grey dashed 
line). 

the spectral features and redshift determination in § 14.11 

2.2. GRB101219A 

GR B 101219A was detected by SwiftlBNY ('Gelbor d et all 
120101) and Konus-Wind (G olene tskii et al. 2010b5 on 
2010 December 19.105 UT. BAT localized the burst at 
a ground-calculated position of RA=04^58™20.7" and 
Dec=-02 °31'37.1" with a 1.0 ' radius uncertainty (90% con- 
tainment; "KrimrnetaDllOlOal). The 7-ray light curve exhibits 
a double-peaked structure with Tgo = 0.6 ± 0.2 s (15-350 
keV) and U = (4.6 ±0.3) x 10"^ erg cm'^ (15-150 keV; 
iKrimm et alj|2010al) . Konus-Wind observations determined 
£peak = 490!^g ^ keV and f-y = (3.6 ±0 .5) x lO-*" erg cm'^ 
(20- 10"* keV: iGolenetski i et al. 2010b). Based on the short 
duration and high iipeak, GRB 1012 19A can be classified as a 
short, hard burst. The 7-ray properties are listed in Table [T] 

2.2.1. X-ray Observations 

SwiftfXRT began observing the field at 6t = 40 s 
and detected a fading, uncatalogued X-ray source at 
RA=04''58"'2 0.49^ and Dec=-02° 32^23.0" with fi nal accu- 
racy of 1.7" (iGoad et al J 120071: lEvans et all 12009 1: Table [B- 
We re-bin the XRT data and extract the best-fit spectrum for 



GRB 101219A as described in EH] We utilize the PC data 
set, St = 70-10"^ s, where there is no evidence for spectral 
evolution. We find an average best-fitting spectrum charac- 
terized by r = 1.8 ± 0.1 and A^H.im = 6.6^^1 x lO^'cm'^ at 



z = 0.718 (C-stat^ = 0.97 for 211 d.o.f.; see §4Jfor redshift 
determination) in excess of the Galactic absorption, A^h.mw = 
4.9 X 102°cm-2 (KalberlaetiD 12005). Our best-fit parame- 
ter s are consistent w ith the automatic spectrum fit produced 
by iGelbord & Grupel (l20Tol) . Applying these parameters to 
the XRT data, we calculate the count rate-to-flux conversion 
factors, and hence the unabsorbed fluxes (Figure|3). 

In addition, we obtained a 20 ks observation with the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S) on-board the 
Chandra X-ray Observatory starting at St = 4.1 days. We ana- 
lyze the Chandra data with the C I AO data reduction package. 
In an energy range of 0.5-8 keV, we extract 4 counts in a 
2.5" aperture centered on the XRT position, consistent with 
the average 3a background level calculated from source-free 
regions on the same chip. We take this count rate of < 2 x 10"^ 
counts s"' to be the 3a upper limit on the X-ray afterglow flux 
at St « 4.1 days. Applying the spectrum extracted from the 
XRT data, this count rate corresponds to < 1.9 x 10"'^ erg 

-2 -1 

cm s . 

The X-ray light curve is characterized by a steep decay and 
a short plateau for St < 200 s, followed by a steady decline to 
the end of XRT observations at St « 10"* s. To quantify this 
decay rate, we utilize the s ingle-parameter x^-minimization 
method described in ^2.1.11 Excluding the XRT data at St < 
200 s and the late-time upper limits, the best-fit power law 
index is ax = -1.37 ±0.13 (xl = 1.1 for 5 d.o.f.). The full 
X-ray afterglow light curve, along with the best-fit model is 
shown in Figure [3] 

2.2.2. Optical/NIR Observations and Afterglow Limits 

UVOT commenced observations at St = 67 s. Over St = 
67-5500 s, no corresponding source was found within the 
XRT position to a 3g limit of > 21.4 in the white filter 
(iKuin & Gelbordll20T0l) . 

We observed the field of GRB 1012 19A in both r- and 
/-bands with GMOS on Gemini-South, and in /-band with 
FourStar, starting at St « 0.96 hr (Table|2]i. We detect a single 
extended source within the XRT error circle in all filters. To 
assess any fading, we obtained additional observations in the 
W-bands at St « 0.2 d (Table|2l). Digital image subtraction be- 
tween these epochs does not reveal any residuals (Figure |4), 
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Table 2 

Log of Optical/NIR Afterglow and Host Galaxy Photometry 



GRB 


Date 


St 


Telescope 


Instrument 


Filter 


Exposures 


Qfwhm 


Afterglow" 


K 


Ho,st° 


^A.MW 




(UT) 


(d) 








(s) 


(") 


(AB mag) 




(AB mag) 


(mag) 


D D 1 AAAO ^ A 

LiKd IUUoZjA 


zUiU Jun ZD. Zoo 


0.52 


Gemini-S 


GMOS 


r 


5 X 120 


2.31 


> 22.6 


< 3.3 


ZZ. /□ ± U.ZJ 


0.027 




n Tun "^ni 


U.J J 


Gemini -S 


GMOS 




3 X 120 


1 Q1 


> 22 7 


^ z.v 


oo 1 n -t- n t; 
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22.04 ± 0.07 
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22.59 ± 0.13 


0.027 
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Magellan 




i 


2 X 240 


0.83 






> 22.7 


0.679 


GRB 101219A 


2010 Dec 19.15 


0.04 


Gemini-S 


GMOS 


i 


9 X 180 


0.66 


> 24.9 


< 0.40 


23.20 ± 0.1 1 


0.097 




2010 Dec 19. 16 


0.05 


Magellan 


FourStar 


J 


25 X 60 


0.46 


> 23.6 


< 1.36 


22.43 ± 0. 13 


0.041 




2010 Dec 19.17 


0.07 


Gemini-S 


GMOS 


r 


9 X 180 


0.80 


> 24.9 


< 0.40 


23.83 ± 0.26 


0.131 




2010 Dec 19.20 


0.09 


Gemini-S 


GMOS 


i 


9 X 180 


0.69 


> 24.9 


<0.40 


23.40 ±0.09 


0.097 




2010 Dec 19.27 


0.16 


Gemini-S 


GMOS 


r 


12 X 180 


0.67 


> 25.1 


< 0.34 


23.73 ±0.10 


0.131 




2010 Dec 19.30 


0.20 
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GMOS 


i 


12 X 180 


0.67 






23. 19 ±0.08 


0.097 




2010 Dec 28.16 


9.05 


Gemini-S 


GMOS 


r 


12 X 240 


0.65 






23.95 ±0.05 


0.131 




2011 Jan 12.15 


24.05 


Magellan 


LDSS3 




6 X 180 


0.68 






23.22 ±0.16 


0.072 




2011 Jan 12.17 


24.06 


Magellan 


LDSS3 


g 


5 X 180 


1.05 






24.57 ± 0.08 


0.189 




201 1 Dec 07.24 


353.1 


Magellan 


FourStar 


J 


15 X 60 


0.56 






22.11 ±0.19 


0.041 




201 1 Dec 07.25 
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Magellan 


FourStar 


Ks 
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0.44 






21.55 ±0.21 


0.017 


GRB 110112A 


2011 Jan 12.18 


0.64 


WHT 


ACAM 


i 


2 X 300 


1.10 


22.77 ±0.29 


2.84 ±0.75 




0.104 




2011 Jun 27.83 


166.2 


Magellan 


LDSS3 


i 


5 X 240 


0.94 






> 24.7 


0.104 




2011 Jun 27.83 


166.3 


Magellan 


LDSS3 


r 


3 X 360 


1.11 






> 25.5 


0.140 




2011 Jul 28.46 


197.3 


Gemini-N 


GMOS 


i 


15 X 180 


0.61 






> 26.2 


0.104 



Note. — Limits correspond to a 3(T confidence le vel. 

" These values are corrected for Galactic extinction iSchlaflv & Finkbeinei|2Ql'll) . 

Only applies to approximately half of the eiTor circle. 
^ Magnitudes for SI and S4, respectively. 

SI is blended with a neighboring bright star (Figure fT6l so we cannot perform photometry. 



^Sfia 10121 9A:-r^s{^:yt^ 




Figure 4. Gemini-South/GMOS i-band observations of the host galaxy of GRB 101219A. The XRT error circle has a radius of 1 .7" (90% containment; black). 
An additional I'-band observation at St = 2.2 hr, adds no additional constraints so is not shown here. Left: St = 0.96 lir. Center: St = 4.8 lir. Right: Digital image 
subtraction of the two epochs reveals no afterglow to a 3(t limit of > 24.9 mag. 



allowing us to place limits on the optical afterglow of / > 24.9 
mag and r > 24.9 mag at the time of the first epoch for each 
filter: St w 0.96 and 2.2 hr, respectively (Table |2]l. To assess 
the fading on timescales > 1 day, we obtained a third set of 
observations in the r-band at w 9 d. Image subtraction with 
each of the first and second r-band observations also show no 
evidence for fading (Table |2]i. A second set of /-band obser- 
vations at 6t w 350 d and a clean image subtraction with the 
first epoch allows us to place a limit on the NIR afterglow of 
7 > 23.6 mag at (5f = 1 .7 hr Finally, to complement our early 



optical/NlR observations, we obtained imaging of the puta- 
tive host galaxy in the ^z-bands with LDSS3 starting on 201 1 
January 12.15 UT, and in the /Tj-band with FourStar on 201 1 
December 07 . 24 UT. Our limits for the afterglow and photom- 
etry of the putative host galaxy are summarized in Table |2] 

We obtained spectroscopic observations of the host on 201 1 
January 2.25 UT using GMOS on Gemini-North at a mean 
airmass of 1.2. We obtained a set of 4 x 1800 s exposures 
with the R400 grating and an order-blocking filter, OG515 
in the nod-and-shuffle mode, covering 5860- 10200 A at a 
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Figure 5. Swift/XRJ light curve of GRB 1 101 12A. The data (red points) for 
St > 200 s is best fit with a single power law characterized by qx = -1 ■ 10 it 
0.05 (grey dashed line). 

spectral resolution of « 7 A. We used standard tasks in IRAF 
for data reduction, CuAr arc lamps for wavelength calibra- 
tion, and archival observations of the smooth-spectrum stan- 
dard star BDh-28 4211 for flux calibration. We discuss the 
characteristics of the spectrum and redshift determination in 



2.3. GRB110112A 

S wiftl^P^ detected GRB 1 101 12A on 201 1 Januai-y 12. 175 
UT (IStamatikos et alJ[201 ll) . with a single spike with Tgo = 
0.5 ±0.1 s (15 -350 keV ) and = (3.0 ± 0.9 ) x 10"** 
erg cm-2 (15 - 150 keV; iBarflielmv etal] |20TT1) . The 
BAT ground-calculated position is RA=21*'59™33.6' and 
Dec=+26°2 8^10.6'^ wifli 2 .6' radius uncertainty (90% con- 
tainment; Barth elmv et al.ll20ril) . The 7 -ray properties are 
listed in Table [T] 

2.3.1. X-ray Observations 

XRT commenced observations of the field of 
GRB110112A at St = 76 s and located a fading X-ray 
counterpart with a UVOT-enhanced positional accuracy of 
1 .6" radius dEvans et al.ll20TTt iGoad et al.ll2007t lEvans etal.1 
l2009t TablelB. We extract the XRT light curve and spectrum 
in the manner described in Section 12.1.11 requiring a mini- 
mum S/N = 3 for each bin, and use the Galactic absorption 
in the dire ction of the burst of A^hmw = 5.5 x 10^** cm~^ 
(iKalberla e t al. 2005). The fight curve is cfiaracterized by a 
sfiort plateau for St < 200 s, followed by a steady decline 
(Figure |5]l. Performing x^-minimization, we find the XRT 
light curve for 5t > 200 s is best fit with a single power law 
characterized by index ax = -1.10 ±0.05 (xi =1.0 for 17 
d.o.f.). Our best-fitting spectral parameters over the entire 
data set, where there is no evidence for spectral evolution, are 
r = 2.2 ± 0.2 and an upper limit of A^h int < 1 .6 x 10^' cm"^ 
(3cr at z = 0; C-stat = 0.82 for 156 d.o.f.). 

2.3.2. Optical Afterglow Discovery 

UVOT commenced observations at St = 80 s, and no cor- 
responding source was found within the XRT position to 
a 3(7 limit in the white filter of > 21.3 mag using data 
over St = 44 00-6100 s (uncorrected for Galactic extinction; 
iBreeveld &"Stamatikos 201 1.) . 



GRB,iiaii-aA ' 

WHT/AeAM i-bafrd 
3o3 



Magellan/LDSS3 i-bahd 
2011 Jun 27.83 tlT 



2031 Jan12..18UI-'-.^.., .,.^^1^ ^ . 



5 arcsec 




Figure 6. The optical afterglow of GRB 1 101 12A. The XRT error circle has a 
radius of 1.6" (90% containment; black) and the red cross marks the centroid 
of the optical afterglow, with a Icr uncertainty of 0.14" (afterglow centroid 
+ absolute tie to SDSS) and i = 22.77 ±0.29. Left: WHT/ACAM i-band 
observations at St = 0.64 days. Right: Magellan/LDSS3 i'-band observations 
at St = 166 days. 

We obtained /-band observations with ACAM mounted on 
the 4.2-m WilHam Herschel Telescope (WHT) at St = 15.4 
hr In a total exposure time of 600 s (Table |2]i, we detect 
a single source within the enhanced XRT error circle with 
/ = 22.77 ± 0.29 mag, where the zeropoint has been deter- 
mined using sources in common with the SDSS catalog (Fig- 
ure |6]l. To assess any fading associated with this source or 
within the XRT position, we obtained /-band imaging with 
LDSS3 starting on 201 1 June 27.83 UT and no longer detect 
any source within the error circle to / > 24.7 mag, confirming 
that the source has faded by > 2 mag. Therefore, we consider 
this source to be the optical afterglow of GRB 1101 12 A. 

To determine the position of the afterglow, we perform 
absolute astrometry using 108 point sources in common 
with SDSS and calculate an astrometric tie RMS of 0.11". 
The resulting afterglow position is RA=2l''59™43.85' and 
Dec=+26°27'23.89" (J2000) with a centi'oid uncertainty of 
0.09" determined with Source Extractor, which, together with 
the astrometric tie uncertainty, gives a total positional uncer- 
tainty of 0.14". We note that this source's position i s not con- 
sistent with the R = 19.6 ± 0.3 source claimed by IXin et al.1 
(I2OI lb . Furthermore, we do not detect any source at this po- 
sition in any of our observations. 

To perform a more thorough search for a coincident host 
galaxy, we obtained r-band observations with LDSS3 on 
2011 June 27.83 UT and /-band observations with Gemini- 
North/GMOS on 201 1 July 28.46 UT. In these deeper obser- 
vations, we do not detect any sources within the XRT error 
circle to limits of r > 25.5 mag and / > 26.2 mag (Table|2]). 
We further assess the proba bility of potential host galaxies 
outside the XRT position in ^4.31 

3. AFTERGLOW PROPERTIES 

We utilize the X-ray and optical/NIR observations to con- 
strain the explosion properties and circumburst environ- 
ments of GRBs 100625 A, 1012 19A and 1101 12A. We adopt 
the standard synchrotron model for a relativistic blastwave 
in a constant density mediu m (ISM), as expected for a 
non-m assive star progenitor (iSari et al.lll999i iGranot & Saril 
I2OO2I) . This model provides a mapping from the broad-band 
afterglow flux densities to physical parameters: isotropic- 
equivalent kinetic energy (/iK.iso), circumburst density (no), 
fractions of post-shock energy in radiating electrons (e^) and 
magnetic fields (e^), and the electron power-law distribution 
index, p, with A^(7) cx j''' for 7 > 7inin. Since we have opti- 
cal and X-ray observations for these three bursts, we focus on 
constraining the location of the cooling frequency (I'c) with 
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respect to the X-ray band because it affects the afterglow flux 
dependence on /iK.iso and no- For each burst, we determine 
this by comparing the temporal (ax) and spectral {f3x = 1 - F) 
indices to the closure relation a- 3/3/2: for p > 2, if i/c > vx, 
a-3/3/2 = 0, while for < vx, a-3/3/2 = 1 /2. We also infer 
the extinction, Aj)?^' by a comparison of the optical and X-ray 
data. 

3.1. GRB100625A 

From the X-ray light curve and spectrum of GRB 100625 A, 
we measure a temporal decay index of ax = -1 -45 ± 0.08 and 
a spectral index of /3x = -1 -5 ± 0.2, which gives ax - 3/3x/2 = 
0.79 ± 0.34. This indicates that Vc < vx and therefore p = 
2.7 ±0.2. 

From our derived value of A^H^nt ^ 1-7 x 10^' cm"^, we 
infer Ay°^' < 0.8 mag (3cr) in the rest-frame of the burst us- 
ing the Galactic Nh-^o-Av conversion, A^H.intMv ~ 2.0 x 10^' 
jPredehl & Schm itt 1995; Watson 2011). We can also inves- 
tigate the presence of extinction by comparing the X-ray flux 
and the optical upper limit at St w 0.5 d. If we assume a max- 
imum value of l^c max ~ 

2.4 X 10'^ Hz (1 keV) and extrapolate 
the X-ray flux density of « 9 x 10~''/iJy to the optical band 
using (3 = -{p- l)/2 = -0.85 to obtain the lowest bound on 
the expected afterglow flux in the absence of extinction, we 
estimate /v.opt ~ 0.24 /iJy (/ = 25.4 mag). Given the observed 
limit of /v_opt S 2.9 /iJy (/ > 22.7 mag), this does not conflict 
with this lower bound, the afterglow observations are consis- 
tent with no extinction. 

We can therefore use the X-ray data and optical afterglow 
limits to constrain /sK.iso and hq. Assuming that the X-ray flux 
is from the forward shock, we can directly obtain E^.iso by 
dGranot & S ari 2002) 

p4.7/4 1.7 0.7/4 ,^-3 

where iiK.iso.52 is in units of 10^^ erg, and ee and eg are in 
units of 10"', and we have used z = 0.452. The X-ray flux 
density at St « lO'* s is F^^^x ~ 9.1 x lO'VJy (1 keV), and 
therefore ^K.iso « 1.2 x 10**° erg (e, = = 0.1). At z = 0.452, 
E^^i,o w 4.3 X 10^° erg (20-2000 keV from the Konus-Wind 
fluence), which gives a 7-ray efficiency of -q^ k, 0.8. If we 
instead assume gf = 0.1 and = 0.01, then /sk.Iso « 1.7 x 10^" 
erg, and r]^ k, 0.7. 

For v„, < I'opt < Vc (where v„, is the synchrotron peak fre- 
quency), the optical afterglow brightness depends on a com- 
bination of /iK.iso and «(). Therefore, the r/z-band limits on the 
afterglow translate to an upper limit on the physical parame- 
ters, given by 

p5.7/4 „0.5^1.7 3.7/4 ^ ^ ^ , „_3 

^K,iso,52"0 ^e.-\^B-\ ^-5 X 10 , (2) 

where no is in units of cm"^. Assuming = = 0.1 and 
using /iK.iso = 1.2 x 10^° erg, we obtain n < 1.5 cm"^. If we 
instead assume = 0.1 and = 0.01, then no ^ 40 cm"-^. 
For both scenarios, we obtain i^e <^ 4 x 10'^ Hz (> 0.02 keV), 
consistent with our assumption that Vc < vx- 

3.2. GRB101219A 

From the X-ray Ught curve and spectrum, we measure ax = 
-1.37 ±0.13 and = -0.8 ±0.1, which gives ax-l>l3x/2 = 
0.17 ±0.23, suggesting that Vc > vx- The resulting value of 



p is 2.7 ±0.1. We note that the closure relation is consistent 
with the alternative scenario for > 2(t. 

Since the optical afterglow flux may be subject to an ap- 
preciable amount of extinction, as suggested by the intrin- 
sic absorption in the X-ray spectrum ( ^2.2. H . the most re- 
liable proxy for /iK.iso and «o is the X-ray afterglow flux. 
Using the last XRT data point at (5? w 7 x lO-' s, which has 
Fy}i ~ 0.03 /xJy (1 keV), we infer the following relationship 
between ZiK.iso and no, 

„0.5 1.7 3.7/4 _ , n -, (^-3 /^x 
-'^Kaso^Sa^O ^e-\^B.-\ ~i-5XlU , (^i) 

where we have used z = 0.718. At this redshift, we find 
w 4.8 X 10^' erg (20- 10"* keV using the Konus-Wind 
fluence). Assuming is-y.iso ~ ^K.iso, we infer no w 1.3 x 10"^ 
cm"-' for Eg = eg = 0.1. With these values, w 6 x 10'^ 
Hz (250 keV), consistent with our assumption that I'c > vx. 
We note that this assumption is violated for no > 4 x 10"^ 
cm"-'. If instead we use e^, = 0.1 and = 0.01, then we obtain 
no w 9 X 10"'* cm"3 and jy, w 2 x 10''' Hz (80 keV), which is 
again self-consistent, and find this assumption is violated for 
no > 0. 1 cm"^. Therefore, the X-ray data suggest an explosion 
environment wifli no w 10"^ - 10"^ cm"^ for GRB 1012 19A. 

We investigate the presence of extinction intrinsic to the 
host galaxy by comparing the X-ray and NIR observations, 
since the NIR data provide a stronger constraint than the opti- 
cal band. Since the X-ray and NIR bands lie on the same seg- 
ment of the synchrotron spectrum, the spectral slope is given 
by /3nir-x = Px ^ -0.8. At the time of our first /-band obser- 
vations at (5f « 1 hr, the X-ray flux density is 0.06 /iJy, lead- 
ing to an expected 7-band flux density of /v,y ~ 14.7 /iJy (21 
mag). This is above the limit of our observations, < 1.4/iJy 
(> 23.6 mag), indicat ing that Aj ^, 2.5 m ag. Using a Milky 
Way extinction curve ( ICardelli et al.llI989l) . this indicates that 
^host > 4 2 mag in the rest-frame of the burst. In addi- 
tion, using the Galactic relation between Nh and Ay, this 
implies A^H.int ^ 7.5 x 10^' cm"^, which does not necessar- 
ily violate our inferred value from the X-ray spectrum of 
A^H,int = (6.6 ±2.0) X 10^' cm"2. Therefore, the broad-band 
afterglow spectrum requires an appreciable amount of extinc- 
tion. 

3.3. GRB110112A 

From the X-ray light curve and spectrum, we measure ax = 
-1 . 10 ± 0.05 and = -1 -2 ± 0.2, giving ax - 3/2/3x = 0.70 ± 
0.30 indicating Vc < vx. The resulting value of is 2.1 ±0.1. 

From our derived value of A^e.int ^ 1-6 x 10^' cm"^, we in- 
fer Ay°^' < 0.9 mag in the rest-frame of the burst using the 
Galactic relation. We can measure the cooling frequency by 
comparing the X-ray and optical fluxes at St w 0.64 d. At this 
time, /v,z ~6.6 x 10"-' ^Jy and/v^opt w2.8^Jy. Using p = 2.1 
and the location of the optical and X-ray bands, we then esti- 
mate that Vc w 1.6 X 10''' Hz (w 7 X 10"-' keV) which agrees 
with our assumption that Vc < vx. The cooling frequency is 
dependent on a combination of physical parameters and gives 
the constraint: 

£M2«o'e^;4«5.4, (4) 

where we have assumed a fiducial redshift of z = 0.5, the me- 
dian of the observed short GRB population. We then use the 
X-ray afterglow flux at St ~ 0.64 d to determine ZiK.iso by 
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GRB 100625A z = 0.452 
°-^Caii age = 640 Myr 
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Figure 7. LDSS3 spectrum of the eai'ly-type host galaxy of GRB 100625A, 
binned with a 3-pixel boxcar (black: data; blue: enor spectram). Also shown 
is the best-fit SSP template (red; Bruzual & Chariot 2003) with a stellar pop- 
ulation age of 640 Myr at a redshift of z = 0.452 ± 0.002. Fits are perfonned 
on the unbinned data. The locations of the Balmer absorption lines and Ca II 
H&K are labelled. 

<:£524:ii^^:^« 0.023. (5) 

Our final constraint comes from the optical afterglow bright- 
ness, given by 

rn5.I/4 05 II 3.1/4 „ 

^K.iso.52«o e,,-ieB,-i ~0-01- (6) 

Assuming e^, = 0. 1 and z = 0.5, we obtain the solution /iK.iso ~ 
2.5 X 10^° erg, no « 1-5 cm"^ and es ~ 0.08. At this redshift, 
^7,180 ~ 9.5 X 10"^^ erg (determined from the Swift fluence and 
applying a correction factor of 5 to represent ss 1 - 10"* keV). If 
we consider a high-redshift origin for GRB 1101 12A of z = 2, 
then we infer larger energies of iiKiso ~ 3.6 x 10^' erg and 
E-y.iso « 1 .5 X 10^' erg, a lower value of eg « 0.01, and a lower 
density, no w 0.18 cm"-'. In both cases, rj^ « 0.3. 

4. HOST GALAXY PROPERTIES 

4.1. GRB 100625 A 

The XRT position of GRB 100625 A fully encompasses 
a single galaxy, which we call Gl (Figure |2]i. To assess 
the probability that the burst originated from Gl, we calcu- 
late the probability of chance coincidence, Pce(< 5R), at a 
given angular separation, {5R) and apparent magnitude (w) 
for galaxies withi n 15^ (the field of view of our images) of 
the burst position (iBloom et al.ll2002l;lBerge3l2010 a). For Gl, 
we conservatively assume 6R = 3(Txrt ~ 3.4", and calculate 
Pcc{< SR) ~ 0.04. The remaining bright galaxies in the field 
have substantially higher values of Pcc{< SR) > 0.17, and a 
search for galaxies within 5° of the GRB position using the 
NASA/IPAC Extragalactic Database (NED) yields only ob- 
jects with Pec ^ 0.98. From these probabilistic arguments, we 
consider Gl to be the host galaxy of GRB 100625 A. 

To determine the host galaxy's redshift, we fit the LDSS3 
spectrum over the wavelength range of 5200 - 8000 A with 
simple stellar population (SSP) spectral evolution models at 
fixed ages (t = 0.29,0.64,0.90, 1.4 and 2.5 Gyr) provide d as 
part of the GAL AXE V Ubrary ("Bruzu al & CharlotI 120031) : at 
wavelengths outside this range, the signal-to-noise is too low 
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Figure 8. grizJKs photometiy for the host galax y of GRB 10062 5 A (black 
circles). The best-fit model (red squares and line; lMarastonll200^ is charac- 
terized by T ft! 0.8 Gyr and M* fa 4.6 X lO' . 

to contribute significantly to the fit. We use x^-minimization 
with redshift as the single free parameter, and perform the 
fit on the unbinned data. The resulting best-fit redshift is 
z = 0.452±0.002 (x^ = 1.3 for 1861 degrees of freedom), de- 
termined primarily by the location of the 4000 A break and the 
main absorption features of Ca II H&K, H/3, H7, and W5. The 
shape of the break is best fit by the template with r = 0.64 Gyr 
(Figure |7]l, and poorer fits (x?, ^ 2) are found for SSPs with 
younger or older ages. Due to the strength of the 4000 A 
break, deep absorption features, lack of emission lines, and 
old age, we classify this host as an early-type galaxy. 

We do not find an emission feature corresponding to 
[O II]A3727. Using the error spectrum, we calculate the ex- 
pected integrated flux for a 3fT emission doublet centered at 
A = 3727 A with a width of « 10 A. We find an expected up- 
per limit of f[oii] < 4.3 x 10"'^ erg cm"^ s"', which translates 
to L[on] < 2.2 X 10"^" erg s"' at the redshift of the burst. Us- 
ing the standard relation, SFR = (1 .4 ± 0.4) yr"' L[oii].4i 
( lKennicuttll998h . we derive a 3cr upper limit of SFR< 0.3 M© 
yr"' for the host galaxy. 

We use the grizJKs-hwd photometry to infer th e stellar 
population age and mass of the host galaxy with the M arastorJ 
j2005i) evolutionary stellar population synthesis models, em- 
ploying a Salpeter initial mass function and a red giant branch 
morpholo gy. We fix Ajj?^' = mag as inferred from the absence 
of A^H.int 02. 1.1b . z = 0.452 as inferred from the spectrum, and 
metallicity Z = Zq, and allow the stellar population age (t) 
and stellar mass (M*) to vary. The resulting best-fit model 
is characterized by t w 0.8 Gyr, in good agreement with the 
fit to the spectrum, and « 4.6 x lO^M©. The model and 
broad-band photometry are shown in Figure |8] 

4.2. GRB101219A 

The XRT position of GRB 1012 19A fully encompasses a 
single galaxy (Gl; Figure HJ. We perform the same prob- 
ability of chance coincidence analysis described in |4T| us- 
ing 5R = 3axRT and find P,c(< SR) w 0.06 for Gl, while 
the remaining bright galaxies within 5' of the burst have 
Pcci< SR) > 0.23. Furthermore, a search within 5° of the po- 
sition with NED yields only galaxies with P„ (< SR) w 1. We 
therefore consider Gl to be the most probable host galaxy of 
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Figure 9. GMOS-N spectram of the host galaxy of GRB 101219A, binned 
with a 3-pixel boxcar (black). The spectram is corrected for Galactic extinc- 
tion and Aj*"^' = 2.5 mag . The stellar population model has r = 25 Myr (red; 
IBrazual & Charlolll2003) . The [O 11] A3727 and [O III] A5007 emission fea- 
tures are at a common redshift of z = 0.718. Also labeled are the locations of 
the Balmer lines He and H5, and marginal emission features at H/3 and the 
[O m] doublet From [O 11] A3727 we deduce SFR = 16.0 ±4.6 Mq yr"'. 
)Kennicutl|[T99^ . 

GRB 101219A. 

We examine the host spectrum of GRB 1012 19A to deter- 
mine the redshift and physical characteristics of the stellar 
population. We identify two emission features in the co-added 
spectrum at Aobs = 6401.65 A and Aob.s = 8599.50 A that are 
also present in the individual 2D spectra prior to co-addition. 
If these features correspond to [O II]A3727 and [O III] A5007, 
their locations give a common redshift of z = 0.718. Fur- 
thermore, we do not find a common redshift solution for an 
alternative set of features, so we consider the host galaxy 
to be at z = 0.718. In addition, we note the presence of 
marginal emission features at the expected locations of H(3 
and [O III]A4959; however, these locations are contaminated 
by sky line residuals. Finally, we detect absorption at the lo- 
cations of He and HS (Figure |9]l. 

To determine the age and host extinction, we use stellar 
population spec tral templates with fixed ages of r = 5, 25, 100 
and 290 Myr GBruzual & Chariot! |2003|) to fit the continuum; 
ages outside this range do not fit the overall shape of the 
spectrum. We apply correctio ns for both Galactic extinc- 
tion (Ay = 0.16 mag at z = 0; iSchlaflv & Fin kbeineil l20Tlh 
and Aj*/"" at z = 0.7 18 using a Milky Way extinction curve 
dCardelh et al.lll98"9h . The spectrum is best matched with the 
r = 25 Myr template and Ajj?*^' = 2.5 mag. Since there is some 
degeneracy between age and Ay°^\ imposing an older stellar 
population of r = 100 Myr also provides a reasonable match, 
but requires a smaller amount of extinction of Ajj?^' « 2 mag. 
Older spectral templates predict a large break at 4000 A not 
seen in the spectrum, while younger templates lack the ob- 
served absorption lines. Therefore, a likely range of ages for 
the host galaxy is t « 25 - 100 Myr Given the emission fea- 
tures and relatively young age, we classfiy this galaxy as late- 
type. The de-reddened spectrum for GRB 1012 19A, along 
with the 25 Myr model, is shown in Figure |9] 

From the extinction-corrected flux of [O II]A3727, F[oii] « 
8.5 X 10""' erg cm'^ s"', we find L[oii] « 1.1 x lO^^^ erg 
s"' at the redshift of the burst. Using the standard relation 
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Figure 10. grizJKs-band photometry of the host galaxy of GRB 101 21 9A 
(black cii'cles). The best-fit model (red squares and line; IMarastor]|2005i) is 
characterized by AJ,™' 1.5 mag, t ft: 15-25 Myr and M, 1.4 X IO'Mq. 
The age-A'™'*' contours of Icr (blue), 2cr (cyan), and 3(t (red) solutions are 
shown in the inset. 
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Figure 11. Left: Gemini-N/GMOS !-band observations of the field of 
GRB 110112A on 2011 Jul 28.46 UT. The position of the optical afterglow 
is marked by the red cross. The five galaxies with the lowest probabilities of 
chance coincidence are circled and labeled G1-G5. The galaxy with the low- 
est value of P{< 5R) is Gl, located 4.8" from the optical afterglow position. 



dKennicuttll 19981) . we derive a SFR of 16.0 ± 4.6 M© yr"'. 

We use the same procedure described in ^4.1l to model the 
SED of the host galaxy to infer r and M*. We fix z = 0.718 
as inferred from the spectrum, Z = Zq, and allow r, M^,, and 
Ay°^^ to vary. The resulting best-fit model is characterized by 
A^°'' « 1.5 mag, r « 15-25 Myr, and w 1.4 x 10'' M©, 
which is consistent with the parameters derived from the spec- 
trum and afterglow. The broad-band photometry and best-fit 
stellar population model are shown in Figure [TOl 

4.3. GRB110112A 

For GRB 110112A, we do not detect a source in coinci- 
dence with the optical afterglow position or within the XRT 
error circle to a 3(7 limit of /> 26.2 mag in our GMOS-N 
image (Figure fTTTi. To determine which sources in the field 
are probable hosts, we calculate Pcc(< 5R) for 15 galaxies 
within ^ 3' of the GRB position, the field of view of our 




Figure 12. Probability of cliance coincidence, Pcc(< SR), as a function of 
angular distance from the optical afterglow position of GRB 1101 12A. There 
are nine galaxies in the 2' field with Pcc(< SR) < 0.85. The five galaxies 
with the lowest P(< 5R) are labeled G1-G5. The galaxy Gl has the lowest 
probability of chance coincidence P( < 5R) = 0.43. 

GMOS-N /-band image. These galaxies were selected by dis- 
carding noticeably fainter galaxies with increasing SR since 
these objects will have Pee(< SR) ^ I. We find that 9 of 
these galaxies have Pcci< 5R) < 0.85 (Figure [TZt. The two 
most probable host galaxies, Gl and G4 (Figures[TT]and [12}, 
have Pcc(< 6R) = 0.43 and 0.54, respectively, and offsets of 
(57? = 4.8" and 11.1". In addition, we search for bright galaxies 
within 5° of the GRB position using NED, but all additional 
catalogued galaxies have Pcc(< SR) > 0.98. Given the rela- 
tively high values for Pcc(< SR), we do not find a convincing 
putative host for GRB 1 101 12A. 

It is also plausible that GRB 110112A originated from a 
galaxy fainter than the detection threshold of our observa- 
tions. For instance, a « 27 mag host would require 5R < 2.0" 
while a « 28 mag host would require SR<1 .3", to be a more 
probable host than Gl. However, to be a 27-28 mag galaxy 
convincing enough to make a host association {Pcci< SR) < 
0.05) would require a smaller offset of SR < 0.5". We note 
that the lack of po tential host is in contrast to previous "host- 
less" short GRBs (lBergeiil2010al) . The high inferred density 
due to the bright optical afterglow (§3) is suggestive of a high- 
redshift origin as opposed to a progenitor system that was 
kicked outside of its host galaxy. 

5. STELLAR POPULATION CHARACTERISTICS 

Of the 30 short GRBs with host associations (P^c < 0.05; 
Table |3]l, GRB 100625 A is the fifth short GRB associated 
with a spectroscopically-confirmed early-type host galaxy 
jGehrels et al. 2005; Bergeretal. 2005; Bloom etal. 200^ 
12007; Fong et al.ll201 1). near the median redshift of the short 
GRB population (Figure [T3Tl. In contrast, GRB 101219A is 
associated with a z = 0.718 late-type galaxy that is actively 
star-forming with characteristics similar to the majority of 
the short GRB late-type host population (Berger 2009). Fi- 
nally, GRB 1 101 12A joins a growing number of short GRBs 
with sub-arcsecond positions but no obvious coincident host 
galaxy to deep limits of > 26 mag (Berger 2010a), although 
unlike previous events, the case for a large offset is less clear 

Short GRBs with sub-arcsecond positions and coincident 
hosts have a median projected physical offset of 5 kpc 
dFong etal.. .20 10.) which, in the context of a NS-NS/NS-BH 




z 



Figure 13. Redshift distribution of 26 short GRBs with host associations 
and redshifts, classified by type of the host galaxy, either late-type (blue) or 
early-type (orange). 

progenitor, can be interp reted as the result of natal kicks 
and moderate delay times (iFryer et al.lll999t iBelczynski et alJ 
'2006). At the inferred redshifts of GRBs 100625A and 
101219A, the upper limits on the projected physical offsets 
set by the radii of the X-ray positions are < 10.3 and < 12.3 
kpc respectively, which agree with the observed offset distri- 
bution. Assuming a fiducial redshift of z = 0.5, GRB 110112A 
would be located 29 ± 3 kpc away from the closest and most 
probable host galaxy, but this association is much less defini- 
tive {Pcc(< SR) ~ 0.43) than previous host-less bursts (BergeiJ 
2010a). Imaging with the Hubble Space Telescope may en- 
able the detection of a faint coincident host. These offsets are 
in contrast t o long GRBs which have relatively small offsets 
of sa 1 kpc (iBloom et all 120021: iFruchter et al.ll2006l) . From 
afterglow observations, the inferred densities for these three 
events may span a wide range, «o ^ 10"*- 1 cm"^*, while lo ng 
GRBs have values of no > 0.1 cm"^ (ISoderberg et al.ll2006h . 

The host galaxies of GRBs 100625A and 101219A have 
stellar populations that span the observed distribution of short 
GRB hosts. With t « 25- 100 Myr and log(M„/MQ) « 
9.1, GRB101219A is at the low end of both th e short 
GRB age and mass distributions (iLeibler & Berger! 12010). 
This host also has one of the most vigorous star formation 
rates reported for a short GRB host to date (iBergerl |2009[ 
iPerlev et"an 120121: iBerger et al.l 120121) . and an appreciable 
extinction of Ay"'*' > 2 mag. These characteristics match 
more closely with the median p arameters o bserved for long 
GRB h ost galaxies (IChristensen et al.i2004t;iWainwright et all 
2007b|; ILeibler &Bergeril2010h . However, an independent 
study based on the 7-ray properties report a probabili ty that 
GRB 1012 19A is not a collapsar of 94% (Bromber get al.l 
I2OI2I) . Compared to other early-type hosts, GRB 10 0625 A 
has a similar age (0.6-0.8 Gyr) an d SFR limit (Bloom et al.l 
"2006; Berger 2009; LeiblerABerger 2010; Fong et al. 201 \), 
but its stellar mass, log(M, t /Mc^,) ^ 9.7, is the lowest by an or- 
der of magnitude (Leibler & Bergerll2010l) . 

6. HOST GALAXY DEMOGRAPHICS 

To broadly determine and utilize the short GRB host pop- 
ulation, we expand upon the observations presented here and 
investigate the demographics of the bulk of the Swift short 
GRB sample, quantifying the fractions of events that explode 
in different types of environments. We divide the population 
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Short GRB Host Galaxy Moiphologies 






GRB 


-"go 




Type' 90% XRT uncert.'' 


Pcc« SR) 


References 




(s) 




(arcsec) 






Sub-arcsecond localized 


050709 


0.07 / 130 


0.161 


L 


3 X 10"' 


1-3 


050724A 


3 


0.257 


E 


2 X 10"^ 


4-5 


051221A 


1.4 


0.546 


L 


5 X ur' 


6-7 


060121 


2.0 


< 4.1 


9 


2 X 10"' 


8-9 


060313 


0.7 


< 1.7 


7 


3 X 10"' 


10-11 


061006 


0.4 / 130 


0.4377 


L 


4x 10"* 


12-15 


061201 


0.8 


0.111 


H/L 


■■■ /0.08 


9, 16-17 


070429B 


0.5 


0.9023 


L 


3 X 10"-' 


18-19 


070707 


1.1 


< 3.6 


7 


7 X 10"' 


20-21 


0707 14B 


2.0/64 


0.9224 


L 


5 X 10"' 


19, 22-23 


070724A 


0.4 


0.457 


L 


8 X 10"* 


24-25 


070809 


1.3 


0.473 


H/E 


■■■ /0.03 


9, 26 


071227 


1.8' 


0.381 


L 


0.01 


27-29 


080503 


0.3/170 


< 4.2 


H/? 


■■■ /O.l 


9, 30-31 


080905A 


1.0 


0.1218 


L 


0.01 


32-33 


081226A 


0.4 


< 4.1 


7 


0.01 


34-35 


090305 


0.4 


< 4.1 


H/? 


■■■ /0.06 


9, 36 


090426A 


1.3 


2.609 


L 


1.5 X 10"* 


37-38 


090510 


0.3 


0.903 


L 


8 X 10"' 


39-40 


090515 


0.04 


0.403 


H/E 


■■■ /0.15 


9,41 


091109B 


0.3 


< 4.4 


7 




42-43 


1001 17A 


0.3 


0.915 


E 


7 X 10"^ 


44-45 


110112A 


0.5 


< 5.3 


H/? 


0.43 


46, This work 


111020A' 


0.4 




7 


0.01 


47-48 


lllinA^' 


0.5 


1.3 


L 


0.02 


49-50 


XRTonly 


050509B 


0.04 


0.225 


E 3.8 


5 X 10"' 


51-52 


050813'' 


0.6 


0.72/1.8 


E/? 2.9 




53-57 


051210 


1.3 


> 1.4 


? 1.6 


0.04 


14,58 


060502B 


0.09 


0.287 


E 5.2 


0.03 


59-60 


060801 


0.5 


1.130 


L 1.5 


0.02 


61-62 


061210 


0.2 / 85 


0.4095 


L 3.9 


0.02 


14, 63 


061217 


0.2 


0.827 


L 5.5 


0.24'' 


14,64 


070729* 


0.9 


0.8 


E 2.5 


0.05 


65-66 


080123 


0.4/115 


0.495 


L 1.7 


0.004 


67-68 


100206 A 


0.1 


0.4075 


L 3.3 


0.02 


69-70 


100625A 


0.3 


0.452 


E 1.8 


0.04 


71, This work 


101219A 


0.6 


0.718 


L 1.7 


0.06 


72, This work 



Note. — " Swift 15—150 keV. For bursts with extended emission, both the duration of the prompt spike and the duration including extended emission are reported. 
Upper limits on redshift are based on the detection of the UV/optical afterglow and therefore the lack of suppression blueward of the Lyman limit (Aq = 912 A) or Lyman-a line 
(Ao = I216A). 

' L=late-type, E=early~type, ?=inconclusive type, H="host-less". For each host-less burst, we also list the type of the galaxy with the lowest Pec IBergetl2010al and this work). 

Only listed for XRT bursts. ( Goad et al. 2007; Evans et al. 2009) 

Evidence at the 4cr level for extended emission is reported to St ^ 100 s. 
^ Bursts with no optical afterglow, localized by Chandra. 

* Bursts with galaxy type classifications based on extensiv e broad-band photometr y jLeibler & Bergedl2010t iMargutti et alj|20129) . In particular, the host of GRB 070729 has an 
infened age 0.98 Gyr) and stellar mass 4 X 10'° Mq ; ILeibler& Bergei f2O10) more consistent with an early-type designation. 

There exists disagreement in the literature regarding the association of GRB 050813 with an early-type cluster galaxy 31 z = 0.72 )Bergej|2005tlFolev et aU2005t IProchaska et alj 
12006 ) or a high redshift cluster at = 1 .8 I Berger 2006); thus, we only display this burst for completeness but do not include it in our demographics. 

' Despite the relatively high P^.^ , all suiTounding galaxies have P^c of order unity I Berger et al. 200" ^. 

References: (1) Villasenor et al. 2005; (2) Fox'et al. 2005; (3) Hjoith et al. 2005b; (4) Krimm et al.'2005; (5)'Berger et al."2005!; (6)'Cummings et al."2005!; (7)'Soderberg et al.l200g; 
(8) de Ugarte Postigo et al. 2006; (9) Berger 2010a; (10) Markwardt et al. 2006; (11) Roming et al. 2006; (12) Urata et al. 2006 (13) Schadv et al. 2006; (14) Berger et al. 2007; (15) 
D^Avanzo et al. 2009; (16) Marshall et al. 2006; (17) Stratta et al. 2007; (18) Markwardt et al. 2007; (19) Cenko et al. 2008; (20) Gotz et al. 2007; (21) Piranomonte et al. 2008; (22) 
Kodakaetal. 2007; (23) Racusin et al. 2007; (24) Ziaeepour et al. 2007; (25) Berger etal. 2009; (26) Marshall et al. 2007; (27) Sato et al. 2007b; (28) D'Avanzo et al. 2007; (29) 
Sakamoto et al. 2007; (30) Mao et al. 2008; (31) Perley et al. 2009; (32) Pagani et al. 2008; (31) Rowlinson et al. 2010; (34) Krimm et al. 2008; (35) Nicuesa Guelbenzu et al. 20 jl; 
(36) Krimm et al. 2009; (37) Antonelli et al. 2009; (38) Levesque et al. 2010; (39) Hoversten et al. 2009; (40) McBreen et al. 2010; (41) Barthelmy et al. 2009; (40) Gates et al. 200g 
(43) Levan et al. 2009; (44) de Pasquale et al. 2010; (45) Fong et al. 201 1; (46) Barthelmy et al. 201 1; (47) Sakamoto et al. 201 1; (48) Fong et al. 2012; (49) Sakamoto et al. 2012; (50) 
Igargutti et al. 2012a; (51) Gehrels et al. 2005; (52) Bloom et al. 2006; (53) Sato et al. 2005; (54) Berger 2005; (55) Foley et al. 2005; (56) Berger 2006; (57) Prochaska et al. 200^ 
(58) La Parola et al. 2006; (59) Sato et al. 2006a; (60) Bloom et al. 2007; (61) Sato et al. 2006b; (62) Berger 2009; (63) Cannizzo et al. 2006; (64) Ziaeepour et al. 2006; (65) Sato etju 
[2007a (66) iLeibler & Bergea2010e (67) iUehara et al..2008i; (68) ,Ukwatta et aU2008i ; (69) iKrimm et ali2010bi ; (70) iPerlev et aU201 It (71) .Holland et aL,2010tt (72) j Krimm et'all 
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Figure 14. Distribution of sliort GRB environments, according to Table [3] 
The fractions of late-type (blue), early-type (orange), host-less (green) and 
inconclusive (yellow) environments are shown. Top: The distribution of 25 
shoit GRBs with sub-arcscond locahzation are divided into all four categories 
(left), and the 6 host-less bursts are each assigned to their most probable host 
galaxy (right; Berger 2010a and this work). Middle: Our full sample, includ- 
ing 1 1 short GRBs with XRT localizations and probable hosts, is divided into 
all four categories (left), and with the 6 host-less bursts assigned (right). Bot- 
tom: Distiibution of our sample for which ther e is no evidence for ex tended 
emission (left) and for which Pf^c > 0.9 (right: iBromberg et al.l2012h . 

into four host galaxy categories: late-type, early-type, incon- 
clusive (coincident hosts that are too faint to classify as late- 
or early-type), and "host-less" (lack of coincident hosts to 
> 26 mag). All late- and early-type designations are based 
on spectroscopic classification, with the exception of two 
hosts, GRBs 070729 and 1 1 1 1 17A, which are based on well- 
sampled br oad-band photom etry (Table |3] iLeibler & Bergerl 
H)10; Marg utti et all2012ah . 

We then use our classifications to examine the relative rates 
of short GRBs detected in early- and late-type galaxies. In 
the absence of observational selection effects, if the overall 
short GRB rate tracks stellar mass alone, the relative detec- 
tion rates in early- and late-type galaxies should match the 
distribution of stellar mass, which is roughly equal at Z ^ 
dkochanek et al.. .2001; .Bell et aL .2003: .Driver et ah l2007h 



and shows little evolution to z 1 dllbert et al]|2010l) . On the 
other hand, if the short GRB rate depends on a combination of 
stellar mass and star formation, as in the case of Type la super- 
novae ( Sullivan et al. 2006), we expect a distribution skewed 
toward star-forming galaxies, with a late-to-early-type ratio of 
>1:1. 

6.1. Environment Fractions 

We first analyze the subset of bursts with sub-arcsecond lo- 
calization because they have the most unambiguous associa- 
tions. Of the 68 short GRBs detected with SwiJ^ as of May 
2012, there are 25 such events (Table|3), 2 o f which have been 
locahzed with Chandra (GRB 111 020A: 'Fong et al.' '2012|; 
GRB 1 1 1 1 17A: lMargutti et al.ll2012 a; Sakamoto et al. 2012), 
an alternative route to sub-arcsecond positions in the absence 
of an optical afterglow. This population is divided as follows: 
11 (44%) originate in late-type galaxies, 2 (8%) are in early- 
type galaxies, 6 (24%) have hosts of inconclusive type, and 6 
(24%) are host-less (Berger 2010a and this work; Figure [l4l 
and Table m. From probability of chance coincidence argu- 
ments, we ca n assign the 6 h ost-less GRBs to a most probable 
host galaxy. iBergeil (I2010al) investigated 5 events, finding 2 
which likely originated in early-type hosts (GRBs 070809 and 
090515), 1 with a late-type host (GRB 061201), and 2 with 
hosts of inconclusive type (GRBs 080503 and 090305). We 
have shown that the remaining host-less burst, GRB 1101 12A 
lacks an obvious host galaxy 04.3l l, and we classify it as in- 
conclusive. 

Accounting for these host-less assignments in the distribu- 
tion of galaxy types, we do not find a substantial change in the 
relative fractions (Figure [l4t. Considering the 16 bursts with 
definitive host types, the late-to-early-type ratio is 3:1 which 
deviates from the expected 1 : 1 distribution if the short GRB 
rate depends only on stellar mass. Using binomical statistics, 
we test the null hypothesis of a distribution that is intrinsically 
1 : 1 and find that the observed ratio has a p-value of only 0.04, 
indicating that the null hypothesis is disfavored (Table |4|i. 

Because the optical afterglow brightness d epends on the cir- 
cumburst density, «o dGranot & Sarill2002l) . the requirement 
of an optical afterglow for precise positions (with the excep- 
tion of the two bursts localized by Chandra) may affect the 
relative rates of short GRB detection in early- and late-type 
hosts if there is a correlation between average density and 
galaxy type. To assess this potential effect, we broaden our 
analysis to include bursts with a single probable host galaxy 
(Pcc(< SR) < 0.05) within or on the outskirts of XRT eiTor 
circles. This sample comprises 1 1 additional eventfl with lo- 
calizations of 1.5-5.5" in radius (90% containment; TableO, 
bringing the total sample size to 36 bursts. Since we require 
sub-arcsecond localization for a burst to be classified as host- 
less, the relative fraction of these events is artificially diluted 
by the addition of bursts with XRT positions (Figure [T4li. 

Assigning the host-less bursts to their most probable 
host galaxies, we recover a similar distribution to the sub- 
arcsecond localized sample: w 50% late-type, ss 20% early- 
type, and « 30% inconclusive, (Figure[T4]and Table|4|i. Based 
on the 26 bursts with early- and late-type designations, this 
gives a late-to-early-type ratio of 2.3:1 and a low p-value of 
0.04 for the null hypothesis that this distribution is drawn 

We note that two of the bursts in our sample, GRBs 050709 and 060121, 
were first discovered by the High Energy Transient Explorer 2 (HETE-2) 
satellite. 

We exclude GRB 050813 from our sample; see Table|4] 
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Table 4 

Shoit GRB Environment Distributions 





Sample 


Late-type 


Early-type 


Inconclusive 


Host-less 


Total 


L:E ratio^' 


Pi,i„„„,(>L:E)'' 


Reject 1 : 1 distribution?^ 


Sub-arcsec 




1 1 (44%) 


2 (8%) 


6 (24%) 


6 (24%) 


25 


5.5:1 


0.01 


Yes 


Sub-arcsec 


, Host-less assigned 


12 (48%) 


4(16%) 


9 (36%) 




25 


3:1 


0.04 


Yes, marginal 


Sub-arcsec 


-^XRT 


17 (47%) 


6 (17%) 


7 (19%,) 


6(17%) 


36 


2.8:1 


0.02 


Yes 


Sub-arcsec 


+ XRT, Host-less assigned 


18 (50%) 


8 (22%) 


10 (28%) 




36 


2.3:1 


0.04 


Yes, marginal 


Sub-arcsec 


+ XRT, All Inc. are Early-type 


18 (50%) 


18 (50%.) 






36 


1:1 


0.5 


No 


Sub-arcsec 


+ XRT, EE excluded 


12 (43%) 


7 (25%) 


9 (32%,) 




28 


1.7:1 


0.19 


No 


Sub-arcsec 


+ XRT, Fnc > 0.9 


8 (58%) 


3 (21%,) 


3 (21%o) 




14 


2.7:1 


0.11 


No 



Note. — " Late-to-early-type ratio 

p value for finding greater than or equal to the observed L:E ratio from a 1 : 1 binomial distribution. 
^ Assumes a significance level of 0.05. 



from an intrinsically 1:1 distribution. To directly compare this 
2.3:1 ratio to the 3:1 observed ratio for sub-arcsecond local- 
ized bursts, we compute the probability of obtaining a ratio 
<2.3: 1 from a population with a true ratio of 3: 1 using Monte 
Carlo simulations for the binomial distribution. In 10^ trials, 
we calculate a high probability of 0.82, suggesting that there 
is no bias to the environment fractions when analyzing only 
sub-arcsecond localized bursts. 

Next, we address the remaining population of 32 Swift short 
GRBs excluded from the discussion thus far. The majority, 
80%, are affected by observing constraints that are dependent 
on factors completely decoupled from any intrinsic proper- 
ties of the bursts: 15 had Swift re-pointing constraints (Sun or 
Moon) and thus have only 7-ray positions, 7 have XRT po- 
sitions that are highly contaminated (in the direction of the 
Galactic plane or near a saturated star, e.g. GRB 100702A, 
see Appendix), and 4 have XRT afterglows but so far lack ad- 
equate optical/NlR follow-up to determine the presence of a 
host galaxy; thus, we cannot currently distinguish between a 
faint coincident host and a host-less origin for these 4 bursts. 
The remaining 20% (6 events) have no XRT localization de- 
spite rapid Swift re-pointing {6t < 2 min), but have a low me- 
dian fluence of « 2 x 10"*^ erg cm~^ compared to the rest of 
the population with (/y.sgrb) ~ 2 x 10~^ erg cm"^ (15- 150 
keV; Figure \T5[ . Therefore, the lack of detectable emission 
with XRT may be related to an intrinsically lower energy 
scale. In summary, we do not expect the exclusion of these 
32 bursts to have a substantial effect on the relative morpho- 
logical fractions. 

The low observed early-type fraction is likely attributed to 
one of two possibilities: (1) it is more challenging to identify 
early-type galaxies at higher redshifts, and thus a dispropor- 
tionate fraction of the bursts designated as inconclusive are 
in fact early-type; or (2) short GRBs preferentially occur in 
late-type galaxies due to the intrinsic properties of their pro- 
genitors. 

We explore the former option by investigating the inconclu- 
sive population in more detail. Spectral energy distributions 
of early-type galaxies generally lack strong emission lines, 
and the most prominent features, the 4000 A break and the 
Ca II H&K absorption lines, are redshifted out of the range of 
most optical spectrographs forz > 1.5, making spectroscopic 
identifications particularly difficult at these redshifts. How- 
ever, more effective studies selecting for distant early-type 
field galaxies by their photometric optical/NIR colors detect 
a nearly constant num ber of early-types between z w 1-1.5 
(IStanford et al."2004|), with a typical AB colo r of 1 - 4 mag, 
depending on the choi ce of o ptical/NIR filters ( IStanford et alj 
l2004t ITamura & Ohtal 12004 . Of the 10 inconclusive host 



galaxies, 4 have optical/NIR color information but yield only 
poor constraints of < 3 - 5 mag due to NIR non-detections 
and faint optical magnitudes, and 5 lack reported NIR obser- 
vations. The only inconclusive host galaxy with multi-band 
detections, GRB 060121, has R-H 2.4 mag; however, the 
optical afterglow and objects in the vicinity are comparably 
red, sugges ting a z > 2 origin as an explanation for the red 
host color ( Levan et al.ll2006al) . iiT-band imaging to depths of 
> 23 AB mag might enable progress in deducing what frac- 
tion of the inconclusive population is more likely early-type. 
To set an extreme upper bound on the true early-type fraction, 
if we assume that all inconclusive hosts are early-types, the 
projected early-type fraction is 50% (Table|4]i. 

We now turn to the second option, that short GRBs prefer- 
entially originate from late-type galaxies. While the predicted 
demographics of NS-NS/NS -BH merger popu lations are cur- 
rently not well-constrained (Belcz ynski et al. ''2006). we can 
use the observed short GRB population to assess the impli- 
cations for the progenitors. We expect to find roughly equal 
early- and late-type fractions if stellar mass is the sole pa- 
rameter determining the short GRB rate. However, we only 
observe this for z < 0.4 (6 events; Figure [T3]l. For z > 0.4, the 
late-type fraction is consistently higher, with a late-to-early- 
type ratio of >2:1. These results, along with the previous 
finding that the short GRB rat e per unit stellar mass is 2-5 
times higher in late-type hosts (iLeibler & Bergerll2010b . sug- 
gest that the short GRB rate is dependent upon a combina- 
tion of stellar mass and star formation. In the context of NS- 
NS/NS-BH mergers, if the delay times of the systems which 
give rise to short GRBs are very long (> few Gyr), we would 
expect a dominant population of early-type hosts at z ^ 0. In- 
stead, the current demographics show a preference for late- 
type galaxies. Along \ yith the inferred stellar p opulation ages 
from SED modeling (ILeibler & Bergerll2010 !). this suggests 
moderate delay times of < few Gyr For a delay time distribu- 
tion of the form P{t) cx t", this translates to n < -1. We note 
that this r esult is similar to Type la supernovae which have 
n « -1.1 dMaoz et al.ll2010l |2012|) . and is in contrast to pre- 
vious short GRB results which claimed substantially longer 
average delay times of ^ 4 - 8 Gyr for lognormal lifetime dis- 
tributi ons based on smaller numbe rs of events (iNakar et al.l 
,2006; .Zheng & Ramirez-Ruiz 2007i: iGal- Yam et al.ll2008lK 

In summary, we find that unless all inconclusive hosts are 
early-type, the short GRB host distribution is skewed toward 
late-type galaxies, with the most likely ranges for the early- 
and late-type fractions of w 20-40% and « 60-80%, re- 
spectively, for the entire short GRB population. Furthermore, 
for most cuts on the sample we find that the null hypothesis 
of a 1:1 distribution can be mildly or strongly rejected. 
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Figure 15. Fluence, /-y, (15 — 150 keV) versus duration, Tgo for the sub- 
arcsecond localized + XRT sample of 36 Swift short GRBs. Bursts are classi- 
fied by morphological type (Table|3) as late-type (blue), early-type (orange) 
and inconclusive (yellow). Open symbols denote host-less assignments. The 
median 2 X 10"^ erg cm"~ and TgQ f» 0.4 s are labeled. The majority of 



events have /-y 10 



- 10"^ erg cm ^. 



6.2. Comparison with j-ray Properties 

We next investigate whether there is contamination in our 
sample from collapsars by analyzing trends between mor- 
phological type and 7-ray properties. We find that bursts in 
early- and late-type galaxies span the entire distribution of 
observed Tgo for short GRBs, with a median value of 0.4 s 
(Figure [TSl l. Using a Kolmogorov-Smimov (K-S) test, we 
find that the two populations are consistent with being drawn 
from the same underlying distribution (p = 0.43). The claim 
becomes stronger when we compare the combined early-type 
and inconclusive distribution with the late-type distribution 
(p = 0.94). On the other hand, the corresponding K-S tests for 
the fluence distributions (Figure [T5t yield marginal /^-values 
of 0.05, suggesting that bursts associated with early- and late- 
types may not be drawn from the same underlying distribution 

in./7- ,. 

A recent study by Brombe rg et al.l (12012) used the 7-ray 

properties (Tgo and spectral hardness) to derive a probability 
that each event is not a collapsar (Pnc), excluding 8 bursts 
which have reported evidence for extended emission. Of the 
29 bursts that overlap in our samples, 14 have a high probabil- 
ity of not arising from a collapsar (Pnc > 0.9). If these proba- 
bilities are robust, and there is contamination from collapsars 
in our full sample, we would expect the galaxy type fractions 
for the population with Pnc > 0.9 to differ from the overall 
sample. In particular, by including only high-probability non- 
collapsar events, we would presumably be excluding mostly 
late-type galaxies since all long GRBs/collapsars are found in 
star-forming galaxies. Therefore, one would naively expect 
the late-to-early-type ratio to decrease with respect to the full 
sample. However, we find that the late-to-early-type ratio for 
this sample is 2.7:1 (Table|4l Figure [T4li which is higher than 
the 2.3:1 ratio inferred for the sample of 36 short GRBs. 

However, Pnc values are not reported for bursts with ex- 
tended emission. Thus, for a more direct comparison, we 
evaluate the subset of 28 short GRBs without extended emis- 
sion (Figure [T4] |. and calculate a late-to-early-type ratio of 
1.7:1 (Table Interestingly, all bursts with extended emis- 
sion originate in late-type (or inconclusive) galaxies, with the 
exception of GRB 050724A. Since the ratio for the Pnc > 0.9 



population is more skewed toward late-type galaxies with 
2.7:1, the probability of obtaining a >2.7:1 ratio in 14 events 
from an intrinsically 1:7:1 distribution is moderate, 0.37. This 
not only demonstrates no noticeable contamination to the 
short GRB host type distribution when including bursts with 
reportedly high probabilities of being collapsars, but also calls 
into question the reliability or importance of these probabili- 
ties in assessing the true population of short GRBs. 

7. CONCLUSIONS 

We present broad-band observations of three short GRBs: 
GRB 100625 A associated with an early-type galaxy at z = 
0.452, GRB 101219A associated with an active star-forming 
galaxy at z = 0.718, and GRB 1101 12A which has a sub- 
arcsecond localization from an optical afterglow but no co- 
incident host galaxy to deep optical limits, and no convincing 
putative host within 5° of the burst location. These observa- 
tions showcase the diversity of short GRB environments and 
give direct clues to the nature of the short GRB progenitor: 
the moderate physical offsets and low inferred densities can 
be interpreted as evidence for a compact binary progenitor 

We also undertake the first comprehensive study of host de- 
mographics for the full Swift short GRB population, classify- 
ing bursts by their host galaxy type. We emphasize several 
key conclusions: 

1. The sample of sub-arcsecond localized bursts have a 
host galaxy distribution of « 50% late-type, « 20% 
early-type and w 30% of inconclusive type after as- 
signing host-less bursts. The inclusion of bursts with 
Swift/XRT positions and convincing host associations 
(Pcc(< 3R) < 0.05) does not affect the relative fractions. 

2. The observed late-to-early-type ratio is >2:1, and most 
cuts to the sample demonstrate that an intrinsically 1 : 1 
distribution is improbable. The only way to obtain 
equal fractions with the observed events is by assum- 
ing that all inconclusive hosts are early-type galaxies at 

3. The most likely ranges for the early- and late-type frac- 
tions are w 20-40% and « 60-80%. The prefer- 
ence toward late-type galaxies suggests that both stel- 
lar mass and star formation play roles in determining 
the short GRB rate. Furthermore, in the context of the 
NS-NS/NS-BH mergers, the observed short GRB pop- 
ulation is not dominated by systems with very long de- 
lay times, but instead with typical delay times of < few 
Gyr 

4. There is no clear trend between Tgo and host galaxy 
type, while there may be a relationship between and 
host type. When excluding the population of bursts re- 
ported to be likely collapsars (> 90% probability), the 
late-type fraction increases relative to the overall short 
GRB sample, suggesting that these probabilities are not 
reliable in assessing the true population. 

Looking forward, our study has demonstrated that detailed 
observations of short GRB afterglows and environments hold 
the key to understanding the underlying population of progen- 
itors. In particular, we emphasize the importance of deep NIR 
observations to determine the early-type fraction within the 
inconclusive population of hosts, and Hubble Space Telescope 
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observations of short GRBs which lack coincident host galax- 
ies to ground-based optical limits (w 26 mag). A concerted 
analysis of broad-band short GRB afterglows would comple- 
ment this study by providing constraints on the basic proper- 
ties of the bursts (i.e., energy scale, circumburst density), and 
help to determine whether there are any correlations between 



these basic properties and galactic environment. Finally, con- 
straints on theoretical predictions for the relative fractions of 
early- and late-type galaxies which host NS-NS/NS-BH merg- 
ers and their delay time distributions will enable a direct com- 
parison to the observed short GRB population. 



APPENDIX 
GRB 100628A 

GRB 100628A was detected by Swift/BAT and the Anti-Coincidence System on INTEGRAL on 2010 June 28.345 UT with 
Tgo = 0.036 ± 0.009 s (15 - 350 keV), = 2.5 ± 0.5 x lO"** erg cm'^ (15 - 150 keV), and peak ene rgy Epeak = 74 1 ± 1 1 .4 keV. 
The ground-calculated position is RA=15''03"'46.2^ Dec=-3r39'10.2" with an uncertainty of 2.1' (llmmler et al.ll2010h . 

X-ray Observations 

XRT began observing the field at St = S6s and detected an X-ray source in coincidence with the core of a bright galaxy. The lack 
of fading of th is source confirmed by Chandra/ ACIS-S observations at St = 4.4 days suggests an AGN origin (Immler et al. 20ld| 
lBergeiil2010bl) . Furthermore, we use binomial statistics and a 10-pixel region centered on the source to calculate the probability 
of a chance fluctuation, finding a high probability of 15%. Thus, this source is ruled out as the afterglow of GRB 100628A. A 
second candidate afterglow was reported based on 7 counts over 3.8 ks in the time interval 6t = 92-7200 s, which translates to a 
count rate of O.OOn+jJ jjog^ counts s"' (0.3 - 10 keV;'Immler et al.' 2010). UVOT, which commenced observations at St = 90 s, did 
not detect a coincident source to > 20.2 mag (white filter; Immler et al. 2010). 

We re-analyze the same time interval of XRT data and use the ximage routine in the HEASOFT package to measure the 
significance of the source. In a blind search, we find the source has a significance of 2.3cr. Lat e-time XRT an d Chandra 
observations confirm that the source has faded by a factor of ^ 15 from the claimed initial X-ray flux (lBergerll2010bl) . However, 
we do not include this burst in our sample of short GRBs with XRT positions due to the low significance of the initial source. We 
caution against classifying this burst as XRT-localized in future short GRB samples. 



APPENDIX 
GRB 100702A 

Swift/BAT detected GRB 100702A on 2010 July 02.044 UT with Tgo = 0.16 ±0.03 s (15-350 keV) and = (1.2 ±0.1) x 10"^ 
erg cm"2 (15- 150key) at a ground-calculated position of RA=16*'22™46.4' and Dec=-56°32'57.4" with an uncertainty of 1.4' 
in radius (isiegel et al.ll2010h . 

X-ray Observations 

XRT started observing t he field at St = 94 s and identified a fading X-ray counterpart with a final UVOT-enhanced positional 
accuracy of 2.4" (TableFll lGoad et al.ll2007l:lEvan s et al."2009^. UVOT commenced observations at St = 101 s and no source was 
identified in the white filter to a limit of > 18 mag (Siegel et al. 2010). The XRT light curve is best fit with a broken power law 
with decay indices of ax.i = ~0-86;|;o 24 and ax,2 = -5.04-'^q'1'^, and a br eak tim e at St = 202 s ( Evans et al.ll2009b . 

We extract a spectrum from the XRT data (method described in § 12. 1.1b and utilize the full PC data set, where there is no 
evidence for spectral evolution. Our best-fit model is characterized by F = 2.7 ± 0.3 and (4.4 ± 2.0) x 10^' cm~^ in excess of the 
substantial Galactic value, A^h.mw = 2.8 x 10^' cm"^ dKalberla et a l. 2005^. We note that the burst is in the direction of the Galactic 
Center (b = -4.8°) a nd therefore the unc ertainties on A^h,mw are likely larger than the typical 10%. Our results are consistent with 
the automatic fits by lEvans et al.l (l2009h . 

Optical/NIR Observations and Afterglow Limits 

We obtained 7-band observations of the field of GRB 100702A with PANIC at St =1.3 hr (Figure [TDl. We detect 4 sources 
within or near the outskirts of the XRT error circle (S1-S4 in Figure [T6]l. S2 and S3 have ste llar PSFs, while S 1 and S 4 have 
non-stellar PSFs. Previously reported 7-band observations also confirm that S2 and S3 are stars (iNicuesa Guelbenzu et al.l l2012). 
while S 1 and S4 have not been reported in the literaturePl To assess any fading within the XRT position, we obtained a second 
set of 7-band observations at St = 6.1 hr. Digital image subtraction reveal no residuals to a 3(7 limit of 7 > 23.3 mag (Table|2]i. 
We caution that this limit only applies to 2/3 of the error circle due to contamination from the saturated star, S2 (Figure [TSTl. 
In addition, we obtained /-band observations with IMACS at St = 247.3 days and we do not detect any additional sources in or 
around the XRT error circle (Table 0. 

^'^ Our PANIC observations show that source "C" in i j- ci t^u • • * ^ - j vm^ 

E' -7^ — ii !„„,- . . n , ... J cludins Si. The remaining two sources are outside of the XRT error 

icuesa (juelbenzu et al.l <2ul2l) is actually three blended sources, in- circle 
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Figure 16. Magellan/PANIC /-band observations of the host galaxy of GRB 100702A. The XRT error circle has a radius of 2.4" (90% containment; black). 
Left: 5f = 1.3 hr. Center: &t = 6.1 hr Right: Digital image subtraction of the two epochs reveals no afterglow to a 3cr limit of 7 > 23.3 mag. 

Probabilities of Chance Coincidence 

We calculate Pcc{< SR) for SI and S4 to assess either source as a putative host galaxy for GRB 100702A. Source SI is fully 
inside the XRT error circle while S4 lies on the outskirts of the XRT error circle. We perform PSF photometry for both sources 
(Table|2]l, and calculate their probabilities of chance coincidence: Pcc(< SR) « 0.02 for S 1 and P{< SR) « 0.04 for S2 using the 3(7 
XRT position radius of 4.5". This analysis suggests that either source is a likely host for GRB 100702A, and we cannot currently 
distinguish which is more likely. We also note that the significant contamination makes it difficult to exclude the possibility that 
there is a brighter galaxy within the XRT error circle. Therefore, we do not include GRB 100702A in our sample of bursts with 
XRT localization, and consider this field to have observing constraints which prevent more in-depth analysis. 
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